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Abstract - We show how the phenomenon of factorization in a quantum many body system is 
of collective nature. To this aim we study the quantum discord Q in the one dimensional XY 
model in a transverse field. We analyze the behavior of Q at both the critical point and at the 
non critical factorizing field. The factorization is found to be governed by an exponential scaling 
law for Q. We also address the thermal efi^ects fanning out from the anomalies occurring at zero 
temperature. Close to the quantum phase transition, Q exhibits a finite-temperature crossover 
with universal scaling behavior, while the factorization phenomenon results in a non trivial pattern 
of correlations present at low temperature. 
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Introduction. — The concepts of symmetry and corre- 
lations pervade all the modern many body physics 
A system consisting in a very large number of particles 
can be found in different phases and the Landau- Ginzburg 
paradigm of symmetry breaking characterizes the various 
phases in terms of different symmetries. Different quan- 
tum phases are separated by Quantum Phase Transitions 
(QPTs), which are driven by tuning an external control 
parameter h across a critical value he [2] . 

Nevertheless, in the past twenty years it has been under- 
stood that symmetry cannot explain quite all the phases 
of matter [3] . Indeed, different patterns of correlations can 
define different quantum phases featuring unconventional 
transitions . Examples in many-body physics come from 
studies on high Tc superconductors, as well as intermetal- 
lic compounds (heavy fermions) and fractional quantum 
Hall liquids [3||§^. 

Here we analyze quantum correlations in a many 
body system addressing the quantum discord, beyond the 



generic notion of "correlations in a quantum system" 7Jl8 



Quantum correlations are not all captured by entangle- 
ment, because a non vanishing quantum discord results for 
certain separable (mixed) states |^. This study addresses 
some new features of the quantum phases involved in the 
phenomenon of symmetry breaking. Besides the critical 
behavior of the quantum discord at the quantum phase 



transition, the discord displays dramatic changes also at a 
non critical value of the control parameter hf ^ he, where 
quantum correlations vanish, thus producing a factorized 
classical state [9 ,10 . Such factorization can even occur 
within the symmetry broken phase, and it consists in the 
sudden reshuffling of quantum correlations, leading to a 

We show 
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transition in the entanglement pattern 
that this correlation transition at hf is governed by a new 
class of scaling laws, thus signaling a collective nature of 
the phenomenon, even if it is not associated to any sym- 
metry breaking. We speculate that the factorization can 
be associated to exotic quantum phase transitions that are 
not described by symmetry breaking but by a reorganiza- 
tion of entanglement patterns without symmetry breaking, 
like in topological quantum phase transitions 
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We complete our study by detecting how the quantum 
critical and the factorization point affect the quantum dis- 
cord at low-temperature, thus opening the way towards 
actual observations 14 . 

Quantum Discord in the XY model. — The total amount 
of correlations in a bipartite (mixed) quantum state pab 
is given by the mutual information Iab = S{pa) + 
S{pb) — S{pab)i where S{p) ~ — Tr[/)log2 p] is the von 
Neumann entropy. On the other hand, classical cor- 
relations can be defined in terms of the quantum con- 
ditional entropy: S{pAB\{Bk}) = J^kPkSiPAs)^ where 
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Bk) PAB {I ® Bk) is the state of the com- 
posite system AB, conditioned to the outcome Bk (being 
a set of projectors representing a complete measurement 
of the subsystem B) of the measurement, with probabihty 
Pk = Tr[(/ (g) Bk) PAB {I ^ Bk)]- The amount of classical 
correlations C is obtained by finding the set of measure- 
ment on B that disturbs the least the part A, i.e., by 
maximizing C = max^^^j [S{pa) - S{pAB\{Bk})] (here 
we restrict to projective measurements) (7j|8]. The quan- 
tum discord is given by: Q = I — C. In a pure state, Q 
reduces to entanglement. A mixed state though, may con- 
tain quantum correlations that are not accounted in the 
lack of separability (see Ref. [t] for examples). 

The model we study is the spin-1/2 chain with XY ex- 
change couplings in a transverse field h: 
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where ct" (a = x,y,x) are the Pauli matrices on site j, 
7 € (0, 1] denotes the xy anisotropy, while h is the trans- 
verse magnetic field strength. The Hamiltonian H is diag- 
onalized by means of a Jordan- Wigncr transformation fol- 



lowed by a Bogoliubov rotation in momentum space 15 



In the range of 7 we consider hereafter, the system dis- 
plays a continuous QPT at he — 1 of the Ising universality 
class with critical indices i' — z = 1, f3 — 1/8. Because 
of superselection rules, the region \h\ < he is magneti- 
cally ordered and the global Z2 symmetry is broken in the 
thermodynamic limit with non vanishing order parameter 
Qx = (o"^); elsewhere the system is a paramagnet. Al- 
though the ground state of H is generally entangled, for 



specific values of 7 and h it is completely separable 12 
Besides the trivial cases h — and /i— >-oo, where IV'gs) is 
fully polarized, there is a non trivial line of factorization 

-t- 7^ = 1 where, for (cr^) ^ 0, |V'gs) = Uj IV'i) 
within the findings of [9[[l0i 
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hi + ^'"^ = 1 where, for 

This line corresponds to an 
accidental degeneracy of the Hamiltonian [l6|[l7], while 
the entanglement pattern swaps from parallel to anti- 
parallel, with a logarithmically divergent range of bipartite 
entanglement ("entanglement transition 
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In order to compute the classical correlations Cr and 
the quantum discord Qr of two spins A and B at dis- 
tance r, one needs to access the single-spin and the two- 
spin reduced density matrices pA and pABij) (see, e.g., 
Ref. [19] for an explicit expression of the generic two- 
spin matrix in a system with global phase flip symme- 
try). Hereafter we focus on the symmetry-broken ground 
state and on the thermal states of Eq. ([T]). For Z2- 
symmetric states, the non vanishing entries of pA and 
PABir) can be evaluated analytically in terms oi gz = (ct^) 
and go 



hr) 
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In that case we use a fully 



analytic treatment for the quantum discord, obtaining the 
thermal ground state as the zero temperature limit of such 
class of states 20 . For symmetry-broken states, gxzir) 
and gx also need to be accessed. Since the expression 
of gxz{T) is cumbersome 21 , in the latter case we resort 




Fig. 1: Quantum discord Qr{h) between two spins at distance 
r in the XY model at 7 = 0.7 (main plot and left inset) and 
7 = 1 (right inset), as a function of the field h. Continuous 
lines are for the thermal ground state, while symbols denote the 
symmetry-broken state obtained by adding a small symmetry- 
breaking longitudinal field hx = 10~^ and it was computed 
with DMRG in a chain of L = 400 spins; simulations were per- 
formed by keeping m — 500 states and evaluating correlators 
at the center of the open-bounded chain. For 7 — 0.7 and at 
hf ~ 0.714, in the symmetric state all the curves for differ- 
ent values of r intersect, while after breaking the symmetry 
Qr is rigorously zero. At the critical point Qr is non analytic, 
thus signaling the QPT. In the paramagnetic phase, there is 
no symmetry breaking to affect Qr- 



to the numerical Density Matrix Renormalization Group 



(DMRG) for finite systems with open boundaries 22 

Ground state. — As displayed in Fig. [l] the difference 
between the quantum discord Qr for the thermal ground 
state and for the symmetry-broken state is always finite 
in the ordered phase (the mutual information / does have 
the same behavior). Moreover, quantum correlations are 
typically much smaller deep in the ordered ferromagnetic 
phase h < he, rather than in the paramagnetic one h > 
he. Nonetheless, as we shall see, they play a fundamental 
role to drive the order-disorder transition at the QPT, 
where Qr exhibits a maximum, as well as the correlation 
transition aX hf, where Qr is rigorously zero. 

Let us first focus on the quantum critical point, where 
the QPT is marked by a divergent derivative of the quan- 
tum discord (see also 20 23 24]). Such divergence is 
present at every 7, for the symmetry broken state; on the 
other hand, for the thermal ground state, it is not present 
at 7 = 1. A thorough finite-size scaling analysis is shown 
in Fig. [2] proving that z = i/ = 1. For the thermal ground 
state (in the thermodynamic limit), we found that d^Qr 
diverge logarithmically as dhQr ^ In |/i — /ic|, within the 
Ising universality class. 

At the factorizing field hf, all the correlation measures 
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Fig. 2: Finite-size scaling of dhQi for tiie symmetry-broken 
state in proximity of the critical point he- Displayed data are 
for 7 = 0.7. The first derivative of the quantum discord is a 
function of L~'^{h — hm) only, and satisfies the scaling ansatz 
dhQi ~ L'^ X F[L^" {h — h,n)], where hm is the renormalized 
critical point at finite size L and uu = 0.472. We found a 
universal behavior he — /im ~ j^-i-28±o.03 ^j^j-^ j-gspect to 7. 
Inset: raw data of dhQi as a function of the transverse field. 



are zero in the state with broken symmetry (see symbols 
in Fig. [1]); in particular, we numerically found a depen- 
dence Qr ^ [h — hfY X (y^)'^ close to it. Such behavior 
is consistent with the expression of correlation functions 
close to the factorizing line obtained in Ref. 25 , and here 



appears to incorporate the effect arising from the non 
vanishing spontaneous magnetization. The factorization 
phenomenon can be traced also for the thermal ground 
state |26]: it is the unique value of the field where the 
same quantum correlations are present at any length scale 
(left inset of Fig. [T]). We found a rather peculiar depen- 
dence of Qj. on the system size, converging to the asymp- 
totic value Qr^~^°°'' with an exponential scaling behavior 
(see Fig. [s]) . The picture elucidated here suggests the ex- 
istence of a non trivial mechanism leading to the factor- 
ization of the ground state. In 11 18 , it was shown that 
hf marks the transition between two different patterns of 
entanglement. The factorization is thus a new kind of 
zero-temperature transition of collective nature, not ac- 
companied by a change of symmetry, and with a scaling 
law that is new in the panorama of the cooperative phe- 
nomena in quantum many-body systems. We emphasize, 
though, that this transition does not correspond -in this 
model- to a QPT. The factorization occurs without any 
non analyticity in the ground-state wave function \gs{h)) 
as a function of /i, as it is shown by the ground-state fi- 
delity J^{h) = \{gs{h)\gs{h + 5h)\. This quantity (which 
can detect both symmetry breaking and non-symmetry 
breaking QPTs 
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), is a smooth function at /i/. So 
there is no QPT here. Nevertheless, the phenomenon of 
factorization can accompany a topological QPT 



Fig. 3: Scaling of Q\ close to the factorizing field, for 7 = 0.7: 
we found an exponential convergence to the thermodynamic 
limit, with a universal behavior according to e~°''"{h — ft^^'), 

Q ~ 1 [h^f'^ denotes the effective factorizing field at size L, 

while 5{Qi) = Q^^' — QJ^~^°°']. Due to the extremely fast con- 
vergence to the asymptotic value, already at L ~ 20 differences 
with the thermodynamic limit are comparable with DMRG ac- 
curacy. Inset: raw data of Qi as a function of h. The cyan line 
is for L = 30 so that, up to numerical precision, the system 
behaves at the thermodynamic limit. 



speculate that the scaling laws associated to topological 
QPTs are those associated to factorization or other phe- 
nomena of entanglement reorganization. At the level of 
spectral properties of the system, we interpret this result 
as an effect of certain competition between states belong- 
ing to different parity sectors for finite L 



17 ; as these 



states intersect, the ground-state energy density is diverg- 
ing for all finite L (such divergence, though, vanishes in 
the thermodynamic limit). 

Finite temperature. — In order to check how the ob- 
served phenomena are resilient with respect to thermal 
fluctuations, we analyze the quantum correlations at finite 
temperature. The low-temperature behavior is influenced 
by the proximity to critical and factorizing flelds. Close to 
he, the physics is dictated by the interplay between ther- 
mal and quantum fluctuations of the order parameter. A 
V-shaped diagram in the h — T plane emerges, character- 
ized by the cross-over temperature T^ross = \h — hc\^ fixing 
the energy scale [2] . T ^ Tcroes identifies two semiclassical 
regimes. In the quantum critical region T ^ Tcross, quan- 
tum and thermal effects cannot be resolved; here the crit- 
ical properties dominate the physics of the system, even 
at finite temperature. Close to hf and at small T, the 
bipartite entanglement remains vanishing in a finite non 
linear cone in the h — T plane 12 18 . Thermal states. 



though, are not separable, and entanglement is present 
in a multipartite form 28 . In this regime the bipartite 



13 . We entanglement results to be non monotonous, and a rcen- 
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Fig. 4: Finite-temperature scaling of the quantum discord for 
the thermal state close to the critical point. The logarithmic 
scaling is verified : along the critical line we found dhQi\h^ ~ 
a;ln(r) + k, with x = 0.065 for 7 = 0.7. The scaling function 
F shows a data collapse close to the critical point. Inset: same 
analysis for the Ising case (7 = 1); we found an analogous 
scaling behavior with x = —0.0059. 



trant swap between parallel and antiparallel entanglement 
is observed 
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At any temperature T > 0, the state is Z2-symmetric. 
By inspection of Fig.[T] and since Qr is a continuous func- 
tion of T for finite temperatures, we conclude that Qr is 
discontinuous as the temperature is switched on, in all 
the phase h < he- Such discontinuity is also observed 
in the entanglement, even if in that case it is much less 
pronounced and occurs only for h < hf 19 . We now ana- 
lyze how criticality and factorization modify the fabric of 
purely quantum correlations in the h — T plane. 

The first aspect we consider is the thermal scaling close 
to the QPT. The scaling ansatz dnQr = T"" F{T /T^ross), 
where Tcross = \h — hc\'^^ 2|, is verified in Fig. |4] for 
r — 1. Remarkably, the discord scales also in the Ising 
case 7 = 1, despite dhQ is not diverging at T = [2^. We 
then discuss the interplay between classical and quantum 
correlations. In Fig. [5^ we show driQi/Ci] in the h — T 
plane, namely the sensitivity of the relative strength be- 
tween quantum and classical correlations to thermal fluc- 
tuations. We found a 1^-shaped diagram, where the ratio is 
constant along the critical line ft = 1 in the quantum criti- 
cal region T ^ |/i^l|, while it explores the largest changes 
along the crossover region. We remark the asymmetry of 
Fig. [5^ between A < 0, A > 0, taking into account that 
the mechanism leading to the two corresponding semiclas- 
sical regimes traces back to quantum (A > 0) or thermal 
(A < 0) fluctuations [2j. 

We now move to hf, where, for the thermal ground 
state, factorization is marked by the fixed point in Qr (see 
left inset of Fig. [I]). This originates a non trivial pattern 
of correlations: Qr{T) ~ Qr'(T) for any r,r'. We quan- 
tify this behavior by analyzing the average displacement 
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Fig. 5: a) Density plot in the h — T plane of drlQi/Ci] 
close to he', along the critical line the ratio Qi/Ci is con- 
stant with respect to the temperature. The solid straight 
line (T — \h — hc\) marks the boundary of the quantum 
critical region, b) Average quantum discord displacement: 
AQ7 = '^J2Tj^i \Qr.{T)-Qrj (T)|/m(m-l) for m = 5 fanning 
out from the factorizing point ft/ ~ 0.714, where all correla- 
tions coincide at any length scale r, as evidenced Fig. [T] left 
inset. Dashed line is for guiding eyes. Here 7 = 0.7. 



between different Qr fanning out from the fixed point in 
the thermal ground state at h = hf (see AQr in Fig. [5|d). 

Outlook and perspectives. — We studied purely quantum 
correlations quantified by the quantum discord Qr in the 
quantum phases involved in a symmetry-breaking QPT. 
Even if Qr results relatively small in the symmetry-broken 
state as compared to the thermal ground state, it underlies 
key features in driving both the order-disorder transition 
across the QPT at he, and the correlation transition across 
the factorizing field h / . The critical point is characterized 
by a non analyticity of Qr found in the Ising universality 
class. Close to hf, Qr displays uniquely non trivial prop- 
erties: in the thermal ground state quantum correlations 
are identical at all scales; for the symmetry-broken state 
we identified a novel exponential scaling, specific for the 
factorization phenomenon emerging as a new kind of col- 
lective phenomenon occurring in the ground state of the 
system. We remark that this can occur without changing 
the symmetry of the system, as a signature of the fact that 
quantum phases and entanglement are more subtle than 
what the symmetry-breaking paradigm says. Although in 
model Eq.([T]) the factorization happens deep in the sym- 
metry broken phase, its behavior is also particularly rel- 
evant in the context of QPTs involving topologically or- 
dered phases 29 , which are believed to occur because of a 
change of the global pattern of entanglement [l3] instead 
of symmetry. 



At finite temperatures a discontinuity of Qr with T 
is evidenced in the whole ordered phase h < he- We 
expect such discontinuity to be present also for models 
with finite Tc- We proved that Qr displays universal fea- 
tures, and it exhibits a crossover behavior; in particular 
the quantum critical region is identified by the condition 
Q{T)/C{T) = Q(0)/C(0) along the critical line. We have 
found that a non trivial pattern of quantum correlations 
fans out from the factorization of the ground state (where 
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AQr ~ 0), opening the way to experimental detection of 
the phenomenon. 
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